The behavior of f(R) gravity in the solar system, galaxies and clusters 
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For cosmologically interesting f(R) gravity models, we derive the complete set of the linearized 
field equations in the Newtonian gauge, under environments of the solar system, galaxies and clusters 
respectively. Based on these equations, we confirmed previous 7 = 1/2 solution in the solar system. 
However, f(R) gravity models can be strongly environment-dependent and the high density (compar- 
ing to the cosmological mean) solar system environment can excite a viable 7 = 1 solution for some 
f(R) gravity models. Although for f(R) oc —1/R, it is not the case; for f(R) oc -exp(-i?/A 2 #o 2 ), 
such 7 = 1 solution does exist. This solution is virtually indistinguishable from that in general 
relativity (GR) and the value of the associated curvature approaches the GR limit, which is much 
higher than value in the 7 = 1/2 solution. We show that for some forms of f(R) gravity, this 
solution is physically stable in the solar system and can smoothly connect to the surface of the Sun. 

The derived field equations can be applied directly to gravitational lensing of galaxies and clusters. 
We find that, despite significant difference in the environments of galaxies and clusters comparing 
to that of the solar system, gravitational lensing of galaxies and clusters can be virtually identical 
to that in GR, for some forms of f(R) gravity. Fortunately, galaxy rotation curve and intra-cluster 
gas pressure profile may contain valuable information to distinguish these f(R) gravity models from 
GR. 

PACS numbers: 04.50.+h,04.25.Nx,95.36.+x,98.62.Sb 



INTRODUCTION 

The standard theory of gravity (the general relativ- 
ity (GR)) combined with the standard model of particle 
physics failed to explain a wide range of independent ob- 
servations, from the expansion of the Universe, the cos- 
mic microwave background, the large scale structure of 
the universe to galaxy and cluster dynamics. To reconcile 
observations, dark matter and dark energy, as modifica- 
tions to particle physics, were proposed and work surpris- 
ingly well However, equally reasonable in logic, one 
can modify gravity instead to reconcile observations. It 
has been shown that the modified Newtonian dynamics 
(MOND) and its relativistic version Tensor- Vector-Scalar 
theory [2] can replace dark matter at galaxy scales to re- 
produce galaxy rotation curves, and 5-D DGP gravity 
[H and f(R) gravity 0, Q can replace dark energy to 
reproduce the accelerated expansion of the universe. 

Like dark matter and dark energy, viable modifications 
in gravity must pass all sorts of tests from the large scale 
structure of the universe d, 0, H, EJ 0, 11 1 to galaxy 
and cluster dynamics to the solar system tests (SST). 
Unlike GR, which involves metric derivatives no higher 
than second order, f(R) gravity involves also third and 
fourth order derivatives, which caused complications in 
the calculation 

An outstanding question is whether f(R) gravity is 
consistent with SST, which have put stringent constraints 
on the PPN parameter 7 = 1 ± O(10" 4 ) p]|. In this pa- 



rameterization, the Schwarzschild metric takes the form 
of ds 2 = (l-2GM/r+---)dt 2 -(l + 2~fGM/r + ---)dr 2 + 
r 2 dfl. For GR, 7 = 1. Various authors have discussed 
the conditions for f(R) gravity or its extensions to pass 
SST 14|. Utilizing the equivalence between f(R) gravity 



and a special class of scalar-tensor theory, Chiba con- 
cluded that the solution of f(R) — —u 4 /R gravity has 
7 = 1/2 and is thus ruled out by SST [l5(. The equiva- 
lence between f(R) gravity and special scalar-tensor the- 
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ory evoked some controversies (e.g. [16|, |17|). Contro- 
versies also exist in approaches without resort to scalar- 
tensor theory. There are different approaches to vary the 
f(R) action to obtain the field equations. One is the met- 
ric formalism, in which the only independent variable is 
the metric. Another is the Palatini form, in which the 
connection is also an independent variable. Based on the 
Palatini form, [H} concluded that f(R) gravity can in 
general be perfectly consistent with SST. Same conclu- 
sion is reached in the metric formalism by several groups 
(e.g. [l9|, lipj. [2llp. However, based also on the metric 
formalism, [23] confirmed the 7=1/2 solution for the 
f(R) = —[i jR gravity. Furthermore, they pointed out 
that the 7=1 vacuum solution with constant curvature 
can not connect smoothly to the surface of the Sun. This 
work has been extended to a wide range of f(R) gravity 
models [2lj |. 

To clarify this crucial issue, we derive the complete set 
of linearized field equations of the two Newtonian poten- 
tials <f> and ip, for cosmologically interesting f(R) gravity 
models, under the metric formalism. We draw the at- 
tention that results in the Palatini form are in general 
different. The field equations turn out to take simple 
forms under the environments of the solar system, galax- 
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ies or clusters. Based on these equations, we confirmed 
the 7 = 1/2 solution. However, we also find that these 
equations accept the 7 = 1 solution for some forms of 
f(R) gravity, due to the presence of non-negligible inter- 
planetary dust and planets. This large matter density 
in the solar system, when compared to the cosmological 
mean, can significantly suppress the corrections induced 
to GR by f(R) gravity and reduces the field equations to 
the GR limit. Wc further checked on the stability of the 
7 = 1 solution and discussed ways to make it stable. 

Galaxies and clusters are excellent laboratories to test 
f(R) gravity through gravitational lensing, cluster X-ray 
and SZ flux and galaxy rotation curve. We perturb the 
FRW background and derive the equations applicable to 
galaxies and clusters. The derived equations can be ap- 
plied directly to address the above observational prop- 
erties. For some f(R) gravity models, the gravitational 
lensing effect is virtually identical to that in GR. How- 
ever, galaxy rotation curve, cluster pressure profile, the 
relation between the cluster mass, X-ray temperature, 
X-ray luminosity and the SZ flux, are modified, even for 
these f(R) gravity models. These observations then have 
discriminating power for a wide range of f(R) gravity 
models. 



LINEARIZED FIELD EQUATIONS OF f(R) 
GRAVITY 



The f(R) gravity takes the action 

L I (R+ f{R))yf=gd A x , 



(1) 



and the field equation 

(R+f)+fnRv 



~9uv a fR — fR;i 



8nGT u 



(2) 

Throughout the paper, we assume that takes the 
form of ideal fluid with negligible pressure. For f(R) 
gravity models of cosmological interest, the f(R) term 
must vanish at high redshifts and R should approach its 
GR limit, in order not to conflict with early time physics 
such as BBN and CMB. At these high redshift the cos- 
mological density is comparable to the solar system local 
density where SST were carried out. So we would expect 
that the correction induced by f(R) can be significantly 
suppressed. However, this is only true if the GR limit, 
namely R — > 8irGp, is reached. In this paper, we will 
explicitly investigate on the feasibility of such solution. 

In the solar system, galaxies and clusters, we expect 
that the gravitational field is weak and the time varia- 
tion of the field is negligible. Due to different environ- 
ments and different boundary conditions, the linearized 
field equations in the solar system differ slightly from 
that in galaxies and clusters. Thus we will treat the two 
cases separately. 



Field equations applicable to the solar system 

In the solar system, we choose a static metric with the 
proper time 

ds 2 = —g^dx^dx" 

= (1 + 2tjj{x))dt 2 - (1 + 20(x)) dx l > 2 . (3) 

i 

This is just the widely adopted Newtonian gauge in cos- 
mology when dropping the time dependence, where <j> 
and ip are two Newtonian potentials. Since |</>|, -C 1, 
non-vanishing Ricci tensor components are Roo ~ V 2 ip 
and Rij ~ — V </><Jjj- — {<f> + ip) t ij- The curvature scalar 
R ~ -2V 2 -0 - 4V 2 </>. In Eq. H it is safe to neglect 
terms (R + /)</>, (R + f)ip with respect to R + / and 
neglect terms 4>Of R , ipOf R with respect to since 
<fr, ip & r e small. Also, one can approximate the covariant 
derivative fR- t i- t j as the ordinary derivative fjijj, since 

|r^/Ri ff //R,«V~ H < L We then obtain 

Roo(l + f R ) + \(R + f) - Of R = SvrGp , (4) 



Ra(i + fR)-^(R + f)-of R 



(/*),« = , (5) 



+ Ir) - = when i ^ j 



(6) 



Before proceeding to the final results, we point out a 
generic constraint exerted by Eq. [6l for constant curva- 
ture solutions. For this kind of solutions, fR-i-j = fR,ij = 
and {<l>+ip),ij = 0. Thus the coefficient of the r _1 term 
in (j) and tp must be equal (with opposite sign) . In another 
word, the constant curvature solution must have 7=1. 

Since \cf) + ip\ <§; 1, Eq. [5] can be integrated to give 



' + V0(1 + Ir) = ~fR + C r 2 + const. 



(7) 



This result is straightforward to check. From Eq. [7J we 
get (4> + (! + /«) + fa + l/0/fl,« + 0+V>),i/fl,j + (0 + 
VOj'/fl.j = ~fR,ij ■ Since |</>+'0| <SC 1, the last three terms 
in the left hand side is negligible comparing to the right 
hand side, we then obtain (0 + ip) t ij(l + f R ) ~ — /ij,^, 
This is the Eq. [S] to begin with. The integral would 
produce a term axi + bxj in the right hand side of Eq. [7] 
However there is no special direction in the Universe, so 
it vanishes. The term Cor 2 is necessary. It reflects the 
fact that, the flat Minkowski space-time is no longer the 
true background in f(R) gravity. 

Combining Eq. dj El and [3 we obtain 



V 2 (<^-V) 



inGp + 2C 
1 + fR 



(8) 



Eq. [7] and [5] completely determine the gravitational 
field of f(R) gravity, up to a constant Co- Co can 
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be determined by either Eq. [4] [5j or the trace of Eq. 
H (f R - l)R ~2f + 3Df R = SirGp. For example, 
when / = A is the cosmological constant, one can show 
C = //8. 

These equations do not require the condition of spher- 
ical symmetry and can be applied to various environ- 
ments including star forming regions and star or black 
hole accretion disk. However, to clarify the issue whether 
f(R) gravity is consistent with SST, we apply them to 
an idealized case, where a point source with mass M (the 
Sun) is embedded in a background with density pback 
(p = M5d{y) + Pback) to sufficiently large radii. 



The 7=1/2 solution found in the literature [23j, \24\ 
corresponds to a small perturbation induced by the cen- 
tral star (Sun) on top of the homogeneous and isotropic 
vacuum background with curvature scalar Rq. For this 
solution, the perturbation in the curvature scalar, is 
Ri ~ ~2GM/3f RR (Ro)r, where Rq is the curvature 
scalar of the vacuum background. The overall curva- 
ture scalar R = R + Ri <C 87rGpback- Following the 
notation of (30j . we call it the low curvature solution. 
This 7 = 1/2 solution is not only viable for the vacuum 
(Pback = 0, [H, Q) but also hold for the realistic config- 
uration of pback 7^ around the Sun, as explicitly shown 
by [HI]. One can also check that the 7 = 1/2 solution is 
indeed a solution of Eq. [7] and [8] 

A crucial step to obtain the 7 = 1/2 solution starts 
with the trace equation (f R — 1)R— 2f+30f R = —8irGp. 
It can be rewritten in the form of Of R = — 8irGp/3+[2f— 
(fn ~ l)i?]/3. In the limit that 8nGp > \2f - (f R - 
1)R\, one has Uf R ~ V 2 f R = d 2 f R /dr 2 + (2/r)df R /dr = 
—8irGp/3. This is the key equation to reach the 7=1/2 
solution. We call the condition 8nGp ^> \2f — (f R — l)R\ 
as the low curvature condition. As explicitly shown in 
[23L I24I [25I ] , the low curvature condition is satisfied for 
the 7 = 1/2 solution. 

However, the trace equation can have another branch 
of solution. It can be rewritten in the form of R = 8nGp+ 
(f R R-2f+3nf R ). In the limit that \f R R-2f+3Uf R \ « 
8nGp, R ~ 8irGp. We call it the high curvature solution 
and call the condition \f R R — 2f + 30f R \ <C 8nGp as the 
high curvature condition. We will show that, given the 
fact that the local density pback is much higher than p c 
(Pback/Pc 10 6 -10 8 , [HI), the high curvature condition 
can be satisfied and a viable 7 = 1 solution is excited, for 
some forms of f(R) gravity. The low curvature condition 
differ significantly from the high curvature condition [3l| . 
Hence the existence of the 7 = 1/2 solution does not 
invalidate the existence of the 7 = 1 solution. Contrast 
to the 7 = 1 vacuum solution, the new 7 = 1 solution can 
connect smoothly to the surface of the Sun and can pass 
all SST. In the solar environment, it virtually reduces 
to GR, with R ~ 8irGp. For this solution, / — > and 

f R « 1- 



The solution of Eq. [8] under the spherical symmetry is 
1 2GM 



^ = 



1 + f R (r = 0) r 



(9) 



dr 



^Gpback + 2C 2 

-r ar 



1 + /R 



Since now \f R \ <ti 1 we obtain 
2GM 



dr 



([87rGp back + 2C ]r 2 dr) 



Combining with Eq. we obtain 



GM 
r 

GM 



dv 

— I Ai:Gp hack r z dr 



C r 2 



+ const. 



/ dr f A n 2, 2C r 2 
+ / ~2 j 4nGp hack r 2 dr + — — 



const. , 
(10) 



where the constant Co is given by the trace of Eq. O By 
variable transform r —* r = r(l + </>), one can express 
the solution (Eq. [TU|) in the form of the Schwarzschild 
metric and then verifies the 7=1 result. 

Is this solution consistent with all approximations we 
made? To be specific, we discuss two forms of f(ft) 
gravity discussed in the literature, f(R) = —p 4 /R 4} 
and f(R) = -\ 1 H 2 exp(-R/\ 2 H 2 ) Q. For f[R) = 
— /i 4 /f? to drive the late time acceleration, p ~ Ho 
where Ho is the present day Hubble constant. For 
f(R) = -Ai^exp(-i?/A 2 ^) proposed in [H, Ai ~ 1 
and A2 ~ 10 3 can produce virtually degenerate expansion 
rate to that of ACDM cosmology. 

We first check the high curvature condition \f R R — 
2f + 3Df R \ < 8nGp. Since now R ~ 8nGp, \f R \ < 
1, the condition reduces to \Of R \ <^ R ~ 8nGp. The 
high curvature condition requires \d 3 f /dR 3 R/r 2 \ <C 1 
and \d 2 f j dR 2 j r 2 \ <C 1 in general. Here, we have adopted 
the approximations R ~ R/r and R ~ R/r 2 where 

For f(R) = -p A /R, \d 2 f/dR 2 /r 2 \ 
\[d 3 f/dR 3 ][R/r 2 ]\ ~ p*/R 3 r 2 ~ [p c / p) 3 / (H 2 r 2 / c 2 ) ~ 
10 12 (AU/r) 2 > 1. Clearly, f(R) = -p 4 /R model does 
no accept the 7=1 solution and contradicts with SST. 

For f(R) = -\ 1 H 2 exp(-R/\ 2 H 2 ) and A 2 = 
10 3 , exp(-i?/A 2 H 2 ) - 10- 400 . So \d 2 f/dR 2 /r 2 \ ~ 
exp(-i?/A 2 fl" 2 )/(iJ 2 r 2 /c 2 ) - 10~ 370 (AU/r) 2 < 1. One 
can further check \[d 3 f / dR 3 ][R/r 2 ]\ < 1. So for f(R) = 
— ^iHq exp(— R/X2H 2 ). the high curvature condition is 
satisfied and 7 = 1 high curvature solution is accepted. 
Furthermore, the weak field condition is also satisfied 
since \(j)\, \ip\ ~ 10" 8 [AU/r] < 1. The condition \f R \ < 1 
is perfectly satisfied too. Hence-after, we only discuss 
this form of f(R) gravity, unless explicitly specified. 

Is it consistent with SST? Yes. Since Co *C Gp, 
Eq. [TO] is virtually identical to that in GR and causes 
no conflict with SST. Furthermore, in the solar system, 



4 



isfied. 



10- 8 (M/M Sun )(AU/r) » is generally sat- 



|Co r /4>\ ^ 1 is a l so satisfied. We then have 
4> + i) ~ and V 2 (c/>- ?/>) ~ -87rGp. Thus the field equa- 
tions are virtually identical to that in GR and we should 
sense no difference between f(R) gravity and GR. 

This solution relies on the condition that pback ^> Pc- 
Clearly, at sufficiently large distance, this condition is 
violated. However, from Eq. [9] and the definition of 7, it 
is clear that the value of 7 is determined by the behavior 
of f R at r — ► instead of r — * 00. So the 7 = 1 conclusion 
is not affected. 

Is it stable against small perturbations? The stability 
of solutions in f(R) gravity models has been widely dis- 
cussed 



iaHSlHSl- For example, ^ found 
that, the high curvature cosmological solutions may not 
be stable while the low curvature solutions can in general 
be stable. Here we briefly discuss the stability of the two 
solar system solutions and leave the detailed investiga- 
tion elsewhere. 

Perturbing the trace equation (f R — 1)R — 2/ + 
SOf R = SirGp around Rq = 8irGp, we obtain 
Ri + Ri/3fim(Ro) ~ 0. A negative f RR causes ex- 
ponential amplification of the initial perturbation with 
a time scale of t = \3f RR (R )\ 1/2 . For f(R) = 
— XiHq exp(— R/X 2 Hq), f RR (Ro) < and the associated 
t <C 1/Ho, since Rq 3> X 2 H 2 . Thus the high curvature 
solution is highly unstable and thus unphysical for these 
two specific models. The low curvature solution is unsta- 
ble, too. However, the time scale to is now comparable 
to the Hubble time, since i?o ~ Hq. So the problem of 
instability is much less severe. 

However, for f(R) = — XiHq cxp(— R/X 2 Hq), one can 
add another small correction to make it stable in the so- 
lar system. We adopt a form of R 2 , which was originally 
proposed by [26[ to drive inflation. Thus the total f{R) — 
— XiHq exp(— R/ X 2 Hq) + aR 2 /Hq, where a is constant. 
This new form of f(R) must satisfy the following condi- 
tions. (1) We want to preserve the high curvature condi- 
tion. This requires a < {H a r/c) 2 ~ l(T 30 (AU/r) 2 . This 
condition automatically preserve the condition f R (Ro) <C 
1, which requires a <C p c /p ~ 10~ 6 . (2) We require 
fRii(Ro) > to make the high curvature solution sta- 
ble. This requires a > Ai exp(--R§/A 2 i$)/2A§. (3) Fur- 
thermore, we want to preserve the merit of f(R) grav- 
ity to drive the late time acceleration. This is satis- 



fied if Ai J ff 2 exp(-i? FRW /A 2 ff 2 ) > 



aRlwt/Hl This 
requires a <C Ai since the late time FRW background 
curvature scalar i?FRW ~ Hq. (4) We want to preserve 
the success of early universe physics such as BBN and 
CMB under GR and thus we want ai?p RW (zBBN) / <C 
-Rfrw(-Zbbn)- BBN happened at T ~ MeV and thus 
z - 10 10 , where i? F Rw(^BBN) ~ 10 30 i/ 2 . Thus we re- 
quire a < 10~ 30 . Since exp(-i? /A 2 i? 2 ) - lO" 400 , 
all these conditions can be satisfied. Thus as long as 
Ai exp(-i? /A 2 ^)/2A| < a < 10~ 30 , this form of f(R) 
gravity contains a high curvature solution, which is con- 



sistent with early and late time cosmological observations 
of BBN, CMB, H - z relation at z < 2 and SST, while 
being physically stable in the solar system [32I ]. 



Field equations applicable to galaxies and clusters 

Galaxies and clusters are virialized objects embedded 
in the FRW background. Approximately they are static 
in the physical coordinate. So the time scale of the field 
variation is the Hubble time and it is safe to neglect all 
the time derivatives of <fi and tp. For a galaxy or cluster 
at z = l/a — 1, we choose a metric 

ds 2 = (l + 2ijj)dt 2 ~ a 2 (l + 2(/))^2dx 1 ' 2 . (11) 

i 

Eq. [7] is then replaced by [ll| 



i> = -/r(-Rfrw + SR) + /r(-Rfrw) 



(12) 



Here, R = R FRW + SR = 6(d 2 /a 2 + a/a) - 2V 2 ^ - 4V 2 </>. 
Throughout this section, the derivative is with respect to 
the physical coordinate. The term Cor 2 presented in the 
solar system solution vanishes, because now the FRW 
background is the right background. One can simply 
verify Co = by the boundary condition that when r — > 
00, <f> — > 0, tp — > and R — > _Rfrw- 

The Poisson equation (Eq. [H]) is replaced by 



8irG(p m - Pb) 
1 + /r 



(13) 



Here, p m is the matter density of galaxies or clusters and 
Ph is the background matter density. These equations 
have been derived under the quasi-static approximation 
at sub- horizon scales llj . Equations applicable to galax- 
ies and clusters are very similar. A good thing is that 
now the quasi-static approximation is well satisfied at all 
relevant scales in galaxy and cluster environments. 

The gravitational lensing is governed by the combi- 
nation (f> — ip. So, given the same matter distribution, 
the gravitational lensing effect of a galaxy or a cluster in 
f(R) gravity is identical to that in GR, except a change 
in the Newton's constant from G to G/(l + f R ). For 
f(R) to drive late time acceleration, f R > in gen- 
eral, thus the lensing signal will be smaller by a factor 
+ In)- However, for some f(R) gravity models, 
such as f(R) = -Ai^exp(-i?/A 2 ^) + aR 2 /Ho pro- 
posed above, \f R \ <C 1 in galaxy and cluster environ- 
ments, so the difference may not be observable. 

However, galaxy rotation curve and intra-cluster gas 
pressure profile can be significantly different to that in 
GR. The acceleration of a test particle is v = —Vip. The 
gas pressure p is determined by Vp = —pVip. The matter 
density decreases from ~ 10 4 p c close to the center to ~ 
40p c at virial radius. The corresponding variation in f R is 
then comparable to variations in <j>, ip and thus (j> + ip = 
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no longer holds. So in galaxy and cluster environments, 
ijj in f(R) gravity does not follow the Poisson equation, 
as that in GR does. Whether the difference caused is 
observable is currently under investigation. 



SUMMARY 

To investigate f(R) gravity in the solar system, galax- 
ies and clusters, we derive the complete sets of the field 
equations which determine the two Newtonian potentials 
4> and ip, under corresponding environments. In the solar 
system, we found that some f(R) gravity models of cos- 
mological interest do contain physically stable solutions 
which are virtually indistinguishable from GR. 

We predict that gravitational lensing effect of quasi- 
static celestial objects such as galaxy and clusters in f(R) 
gravity is virtually the same as in GR, up to a change in 
the Newton's constant. However, galaxy rotation curve 
and cluster gas pressure profile differ intrinsically from 
that in GR. Thus observations of galaxies and cluster 
dynamics are promising to put useful constraints on /(i?) 
gravity. 
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As seen from later discussion, the high curvature condi- 
tion is equivalent to |D/h| <C 8irGp, while the low curva- 
ture condition is equivalent to — 8nGp. 
The cosmological solution of f(R) = 
-Ai# 2 exp(-i?/Y#o) + aR 2 /Hi with a < 10" 30 
is still unstable. Whether one can find suitable correc- 
tion to f(R) — —XiHq exp(— R/\2Hq) to make both 
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the high curvature solar system solution and the cos- 
mological solution stable requires further investigation. 



